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Abstract 

Bose-Einstein condensation of rotons in helium was considered long 
ago. It was shown that the relative velocity of the normal motion in this 
state must be equal to the Landau critical velocity. We argue that the 
condensation can be attained at a smaller velocity if the temperature is 
low enough. 
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1 Introduction 

As a Bose gas is cooled down, the Bose-Einstein condensation (BEC) emerges to 
avoid a conflict between the statistics and the particle conservation requirement. 
Behaviour of quasi-particles, unlike behaviour of real particles, is not dominated 
by such conflict: their number is not an independent variable, it adjusts itself to 
maximize entropy in equilibrium. Conventional Bose-Einstein condensation of 
quasiparticles is therefore thought to be impossible, an ingenious mechanism [1 
of the condensation for rotons via gap cancellation by a critical (v n — v s = 
v = vl = Ao/Poi where Ao is the roton energy gap and Pq is its momentum) 
superfluid counterflow is requirec£]. 

Actually, even real particles hardly ever conserve exactly. Consider atoms of 
a cold gas in a trap: they can combine into molecules or even evaporate from the 
trap altogether. Nor particle conservation is present in exact relativistic theory. 
The BEC is still possible if it forms much faster than the particle population 
decays. It is therefore necessary to compare all relevant thcrmalization time 
scales. It turns out that any finite counterflow and low enough temperature 
T -c vPo are the conditions favourable for the roton BEC creation. 

Suppose these conditions arc satisfied, and the roton number relaxation is 
much slower than that of their energy and momentum. The roton distribution is 
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According to G.E. Volovik [2J, the possibility of the roton BEC was first investigated 
(unpublished) by him as early as 40 years ago. 
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then characterized by the temperature, velocity, and finite chemical potential fi: 



Arp=(exp^^-l)"\ (1) 

where P is the roton momentum and £ = A + (P — Po) 2 /(2[i ) is its energy. 
Let the z axis run along v direction. The roton distribution argument can be 
expanded in powers of small deviation from its most probable value 
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where f = P v (/i + Po/v), 



9x,y = fx.yV VOV/PO, 
9z = fz- 



(3) 



This expansion is applicable if 

T^hqv 2 . (4) 

Whether the roton BEC state is a Bogolyubov-like gas or a degenerate Bose 
liquid depends on the concentration. Critical concentration (particles per unit 
volume) required for Bose-Einstein condensation is 

(If ^- 

The system is almost ideal (see [3]) if 

N Yli^EL <c i ( 6 ) 

where |4] Vb ~ 10 -38 erg cm 3 is the interaction strength. Combining (J5)) and ([6]) 

we get 



2 5 7r 3 n 6 



T « c(3/2)^ w pr 2 ~ 10 ^- 

This inequality is satisfied as a consequence of ((4|), implying that the conden- 
sation considered is a transition between gaseous phases. 



2 Relaxation 

2.1 Roton number decay 

Suppose initial roton distribution is characterized by some positive chemical 
potential fi. This means that the roton number is greater than that in complete 
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equilibrium. The most important process at low temperature T -C vPq for 
the chemical potential relaxation is the transformation of two rotons into one 
roton and one phonon. Little is known about the transformation probability 
in such collisions. It seems reasonable to assume that the transformation cross 
section (providing the process is allowed at all by the conservation laws) can be 
bounded from above by complete scattering cross section [3] known from the 
experimental viscosity data. 

Momentum conservation for this transformation imposes severe restriction 
on the angle <f> between momenta of the incident rotons. Namely, this restriction 
is reduced to <j) > 2n/3 if the inequality A -C P c is taken into account. Here 
c is the speed of sound. To simplify all assessments below we take T <C A, 
i.e., assume Boltzmann statistics for the rotons: 

N P = exp ^ . (7) 

This is certainly incorrect for the BEC state itself but must provide reasonable 
relative order of magnitude for different relaxation rates. 

As an estimate, not more than about exp(— 2(1 — cos7r/3)wPo/^) fraction 
of all collisions end up with the transformation. For the chemical potential 
relaxation rate this gives (see [J]) 

T » < J-A ex P^^> ( 8 ) 

where N is the total roton concentration. The relaxation rate here is defined 
according to 



/i + T„ V = o. 



2.2 Phonon-roton velocity relaxation 

To find the upper boundary for the phonon-roton relaxation time it is sufficient 
to consider the two-particle scattering of rotons by phonons. Conservation laws 
for this process are (primes denote the finite state) 

P + P = P' + P', 

e + S = e' + £', {) 

where p and e — cp are the phonon momentum and energy. Since P > p, we 
conclude that p w p' and P || P'. 

Let the roton and phonon subsystems be separately in equilibrium. The 
roton and the phonon "bath" velocities will be v and v + 6v respectively. Ve- 
locity relaxation has two distinct time scales T v h,-v\\ an d Tj,ft )t) j_, they correspond 
to <5v || v and to 5v _L v. In linear approximation (Sv <C v) the relaxation is 
described by the equations 

ivy = *«nPpir^J t) ||. 
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where j is the roton momentum density and p r is the roton contribution to the 
normal density p n — dj/dv 

p r± =P N/v. V ' 

The relaxation process is governed by the kinetic equation, in the spatially 
uniform case it gives 



(P-P')(n' + l)n7V P dw 



dp dP 

(2ttH) 3 (2%h) 3 ' 



(11) 



where 



e - p(v + <5v) 
" = 1 exp 1 



the scattering probability is dw — cda and differential cross-section according 
to Ref.g] is 



da 



( PoP 2 
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((n + n', m)(n, n') + — (n, m) 2 (n', m) 2 + a\ dn'. 
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Here n, n', and m are the unit vectors directed along p, p', and P respectively. 
At low temperature most rotons have momentum parallel to the velocity and 
the vector m |j v can be regarded as a constant. The parameter A is given by 



Poc 
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To transform Eq. (|ll|) we employ the usual relation between distribution 
functions 



in' + l)nN P - (n + l)n'N P , « (n + l)n'N P 
and the principle of detailed balance: 



(P' -P,6v) 
T 



ir= 2T * (P ' ~ P)(P ' ~ P ' ^ V)(n + 1)n ' Np ' dw 



dp dP 
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p c 
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For the relaxation rate this gives (lengthy but straightforward transformations 
are omitted) 

_! _ 87r 5 P 2 T 8 / 4 2AP P 2 A 2 \ T 8 P 2 
T P h,v\\ ~ 15c 10 ft 7 p 2 Mo V225 + 15^5 + 35^ + ~3~ J ~ 50 °c 10 /iVW 

x Stt^PqT 8 / 8 2AP P 2 T 8 Po^ 
pfc,«-L l5c 10 ftV V 22 5 45/i c 175^c 2 3 / c 10 fi 7 p 2 ' 1 ' 

2.3 Phonon-roton temperature relaxation 

Let us now find the temperature equilibration rate between the roton and the 
phonon subsystems. The phonon gas temperature is T + 8T and the phonon 
distribution is 

( e-pv N_l 

n= [ exp ¥Tsf- 1 

This function satisfies equality 

£' — £ 

M + l)nNp - (n+ l)n'N P > re (n+ l)n'N P , — ; ST. 
From the conservation laws (J9)) it follows that 

£'-£= P ~ P ° (m,P'-P) = «(m,P'-P). 
Mo 

The energy inflow to the roton subsystem is 

& = / (£ '- £ > ( "' + 1 >»* d »<^<^ = 

^/KP'-m.+iK*-^. 

The phonon-roton temperature relaxation rate defined by 

P r = CrTp^rpST, 

where C r = 3iV/2 is the roton contribution to the specific heat per unit volume, 
can be immediately extracted using obvious similarity between (|15l) and (1121) : 



(15) 
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2.4 Roton-roton velocity relaxation 

Equilibrium within the roton subsystem is reached via the roton-roton collisions. 
Like a two-body problem in classical mechanics these collisions are efficiently 
described (if © holds) in the center of inertia frame. Namely, suppose reduced 
momenta, defined according to ©, of the scattering rotons are gi and g2. After 
a transformation 

G = gi + g 2 , 

g = (gi - g2)/2, 
the net energy of two rotons is given by 

C 2 n 2 

£t + £ 2 - vG - 2fi = 2A + — + 2- (17) 

4mo Mo 

and the conservation laws are simplified to 

G' = G 

9' = 9- 

Accurate definition of the roton-roton equilibration time constants is hardly 
possible (cf. discussion on the establishment of equilibrium of a phonon gas 
in Ref.[5 ]). As an estimate we employ the relations similar to those for the 
phonon-roton relaxation. 



AllV^jjjy ir,-ri,n.) ^P."P,<»" (2Tfi) , (2 , B) 3 
ftxvi - ^ / [P;- Pl , m f Wp ,iV Pl d» J^J^. 



We follow Ref. UJ and adopt the simplest possibility for the roton-roton scatter- 
ing probability 

aw = — 6 [Ex + Ei — t\ + £2 ) 



h v 1 z 1 ' (2tt/i) 3 ' 

For the relaxation time this gives 

_ x 47r|y | 2 e- 2A / T [, , ,2 ( gl + g 2 2 
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and 

_ x _ W*\Vq\*P*ii\I 2 TV* 
Here we substituted the roton concentration according to 

xr = -A/T / p -g?/(2J>o) dF l = 2 g f Mo ~A/T 
J (2tT?j) 3 w(27Tft) 3 

Eventually, using (fTU|). we see that the two rates are equal 



8iV|K.| 2 PoMo /2 r 1/2 



(18) 



2.5 Roton-roton temperature relaxation 

Employing the same approach as in (|16p we express the roton-roton temperature 
relaxation rate through the parallel velocity relaxation rate (|18p 



= v 2 p r{] , = 16iV|Vfr| 2 P oAt ;f p 

r,T (7 T 0^3/2^1/2 ' 



3 Discussion 

As we mentioned above, possibility of the Bose-Einstein condensation of rotons 
depends on the relative magnitude of the different relaxation rates. We begin 
with the notion that within the roton gas the temperature relaxation (TT9|) is 
much faster than the velocity relaxation (|18[) 



r 



r.T 



1 - 2 

> 1. 
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Tr,v 

The latter in turn is always faster than the chemical potential relaxation 



r" 1 

r.v 



T 2vP 



Tfj, 3ir A / 2 y /i u^ T 

This completes the proof that a roton subsystem with an initially narrow 
momentum distribution around some undcrcritical momentum P such that 
< P — Pq < fio^o/Po must pass through a BEC state. 

Whether an initially wide roton distribution can be experimentally con- 
densed by the "phonon cooling" depends not only on the relation between the 
rates of the phonon- roton relaxation ([T3|) . (Q~4|), (fTBj) and of the roton number 
decay ((8|) . The latter has exponential dependence on the temperature and can in 
principle be made arbitrary slow relatively, but at low temperature the phonon- 
roton relaxation is very slow itself and may take too long, therefore demanding 
a very large experimental cell. To overcome this difficulty one could try to con- 
dense rotons at rest in the laboratory frame of reference while the superfluid 
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passes through a capillary. Andreev reflection of rotons [B] at low temperature 
will protect the distribution width in the roton-wall collisions. 

Note that exponentially slow roton number decay (in contrast with the 
power-law temperature dependency of other equilibration rates) is a general 
result and does not depend on exact roton-roton interaction. 

Experimental observation of the roton BEC should be possible by a number 
of techniques: 

• the coherent roton quantum state has finite momentum. Upon BEC 
formation, bulk helium acquires a spatial inhomogeneity [7], a one- 
dimensional density wave. The wavelength is the roton wavelength and the 
modulation direction is the velocity v direction. The periodicity should 
manifest itself as a Bragg peak in the x-ray scattering experiments. Actu- 
ally this roton BEC state is a supersolid as it simultaneously has superfluid 
and crystalline order. Note, that one-dimensional crystalline order is not 
destroyed by Landau-Peierls fluctuations [51 [S] thanks to the true three- 
dimensional superfluid order. 

• another option to probe BEC is to explore excitations of the condensate. 
The roton second sound is well studied in normal roton systems |10j and 
may be used for the condensate detection. 

• the roton distribution can be measured directly by the quantum evapora- 
tion . Delta peak in the distribution function would become an explicit 
proof of the BEC formation. 

Finally, let us remark that the stability analysis performed in Ref.[H] is 
irrelevant for the metastable BEC considered in present paper, because the 
roton number is fixed. 
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